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AN OPTIMALITY CRITERION FOR SIZING MEMBERS OF HEATED 
STRUCTURES WITH TEMPERATURE CONSTRAINTS 


G, Venkateswara Rao,* Charles P* Shore, 
and R. Narayanaswami** 

Langley Research Center 


SUMMARY 

A thermal optimality criterion is presented for sizing members of heated 
structures with multiple temperature constraints. The optimality criterion is 
similar to an existing optimality criterion for design of mechanically loaded 
structures with displacement constraints. Effectiveness of the thermal opti- 
mality criterion is assessed by applying it to one- and two-dimensional thermal 
problems where temperatures can be controlled by varying the material distribu- 
tion in the structure. Results obtained from the optimality criterion agree 
within 2 percent with results from a closed-form solution and with results from 
a mathematical programing technique. The thermal optimality criterion augments 
existing optimality criteria for strength and stiffness related constraints and 
offers the possibility of extension of optimality techniques to sizing structures 
with combined thermal and mechanical loading. 


INTRODUCTION 

Aerospace structures may be subjected to severe mechanical and thermal 
loading. Efficient designs for such structures increasingly require use of 
automated design methods to determine minimum-mass structural arrangements that 
satisfy several mechanical and thermal design requirements. Both mathematical 
programing methods (ref. 1) and optimality criteria (refs. 2 to 4) have been 
proposed for automated design of mechanically loaded structures with strength, 
stiffness, and flutter constraints. 

Recent studies have focused on resizing procedures for structures under 
combined mechanical and thermal loading with strength constraints. In refer- 
ences 5 and 6, steady-state thermal effects were included in an implicit fashion 
through strength requirements; however, temperatures in the structure were not 
updated as a result of structural resizing and significant unknown differences 
can exist between final and initial temperature distributions in the structure. 
Explicit thermal design requirements and updating of temperatures during the 
resizing process were included in a finite-element-based mathematical programing 
method described in reference 7 where membrane element thicknesses and bar ele- 
ment areas were used as design variables to adjust average element temperatures. 
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The present paper presents an optimality criterion which satisfies temper- 
ature requirements and permits updating of structural nodal temperatures during 
resizing. The optimality criterion resizing formulas are similar to those of 
reference 3 for design of mechanically loaded structures with displacement con- 
straints. Only temperatures (not stresses) are considered in the present paper, 
and use of the criterion is illustrated by applying it to several examples with 
both one- and two-dimensional heat-transfer characteristics. An existing 
closed-form solution (ref. 8) is used to assess the effectiveness of the opti- 
mality criterion for one-dimensional problems, and a mathematical programing 
technique described in references 9 and 10 is used for a similar assessment for 
two -d imen si ona 1 problems . 


SYMBOLS 

Values are given in both SI and U.S. Customary Units. The measurements 
and calculations were made in U.S. Customary Units. 

A area of bar element 

A = A/Ao 

a general design variable 

Biot number, h\\^^/v.kQ for bar and hip^yktQ for plate 
b convergence parameter 

c,c-|,C 2 ,C 3 points on bar and plate at which constraints are applied 
conductance matrix 
constraint function 
convective matrix 

convective heat-transfer coefficient 
thermal conductivity 
length of bar 

characteristic length of plate 
length of bar element 
mass of structure 
number of bar elements 


[cj 

F 

[H] 

h 

k 

^b 


I 

m 

Nb 
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Nq heat-generation number, QT-b^/^^QTQ for bar and qip^tgl^ for plaite 

Np number of plate elements 

N-p number of points at which temperatures are constrained 

{q} unit thermal load vector 

q rate of heat input 

{ql'j thermal load vector 

{q}2 convection load vector 

S surface area of plate element 

T temperature 

T = T/Tq 

{t} temperature vector 

{Tq} temperature vector corresponding to {q} 

Tg ambient temperature 

Tg nondimensional ambient temperature 

Tq specified temperature at point c 


Tc»Tci >T q2>^c3 


specified nondimensionalized temperatures at constraint 
points c, c-| , C2, and cj 


t 

V 


x,y 

x,y 

Y 

X 


thickness of plate element 
volume of structure 
= V/AqIi) for bar 
= v/tQip^ for plate 
Cartesian coordinates 

nondimensional Cartesian coordinates, x/ib,p and y/Xfc p> 
respectively 

step-size parameter 

Lagrangian multiplier 

density 
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Subscripts: 
b bar 

i ith bar element 

j jth bar element 

k number of constraint point 

o reference value 

p plate 


OPTIMALITY CRITERION 

An optimality criterion for thermally loaded structures is developed from 
the general treatment presented in reference 3 for problems which can be charac- 
terized by £in objective function m that is dependent on a group of design vari- 
ables a^ 


m = m(ai) 


( 1 ) 


and an equality constraint function that may be implicitly dependent on the 
design variables a^ 


F(ai) - Fc = 0 (2) 


The quantity Fq is the desired value of F(aj^) at point c. 

To obtain the necessary condition for a minimum of m subject to equa- 
tion (2), first append the constraint condition to the objective function with 
a Lagrangian multiplier 


m» = mCa^) + Xj^FCa^) - Fq] 


(3) 


and then equate to zero the first partial derivatives of m* 
the design variables a^ 

3m» 3m ^ 3F 
3aj^ 3aj^ 3a^^ 

Equation (4) is also called the optimality criterion. 


with respect to 


(4) 
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For two-dimensional thermally loaded structures the response quantity F(aj^) 
is the structural temperature and is given by the following partial differential 
equation and the appropriate thermal boundary conditions; 


d 

3x 






+ Q = 


0 


(5) 


where and ky are the material thermal conductivities in the x- and 

y-directions, respectively, and Q is the applied thermal loading. 

Equation (4) is specialized in reference 3 for finite-element represen- 
tations of mechanically loaded structures with deflection constraints. The 
finite-element representation allows both m(aj^) and 3F(aj^)/9aj^ to be 
expressed as sums of individual member contributions. As a result, the simul- 
taneous equations represented by equation (4) uncouple and can be solved by 
simple recursive formulas. 

For structures with both convective and conductive modes of heat transfer 
present, equation (5) can be written in finite-element matrix notation as 


[d](T} = {q}, * {q}2 


( 6 ) 


where [l^ is a heat-transfer matrix and includes terms for both convective 
and conductive heat transfer, {t} is a vector of unknown temperatures, and 
{q}-] and {q}2 applied thermal and convective heat load vectors, respec- 

tively. Since equation (6) is analogous to the corresponding equation for 
displacements in mechanically loaded structures and since use of thermal finite 
elements also permits uncoupling of the simultaneous equations represented by 
equation (4), the optimality criterion for deflection constraints developed in 
reference 3 can be adapted to problems with temperature constraints. Appro- 
priate expressions for the optimality criterion and resizing formulas for bar 
and plate elements are now restated in terms of thermal quantities. Equation 
numbers from reference 3 are also given. For clarity and completeness details 
of the adaptation are presented in the appendix. 


Optimality Criterion and Resizing Formulas for a 
Single Temperature Constraint 

For a single temperature constraint the optimality criterion given by 
equation (4) can be expressed in a form similar to equation (44) of reference 3 
as 


®i 

— = 1 (7) 

Pb 
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for bar elements and 


®j 

— = 1 (8) 
Pp 


for plate elements where 


X{T}T[c]^{TQ}j 

fi/iAi 


and 


®J ■ 


X{T}j[c]j(Tg}j 




(9) 


( 10 ) 


In the expressions for e^^ ®j> X is the unknown Lagrangian multiplier, 

{t} is the vector of nodal temperatures obtained by solution of equation (6) 
for the applied thermal loading, [cj is the elemental conductivity matrix, 

{Tq} is the temperature vector obtained by solution of equation (6) for a unit 
thermal load applied at the constraint point c, and Aj^ are the length 

and area of the ith bar element, respectively, and Sj and tj are the surface 
area and thickness of the jth plate element, respectively. The subscripts i 
and j on the quantities {t}, [C] , and {Tq} indicate that these quantities 

are associated with the ith bar element or the jth plate element. 

The resizing formula given by equation (48) of reference 3 is rewritten 
as 


^i,n+1 



( 11 ) 


for bar elements and 


^J,n+1 



( 12 ) 


for plate elements, where n is the iteration number and y is a paramter that 
controls the step size. 
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Optimality Criterion and Resizing Formulas for 

Multiple Temperature Constraints 

If the structure has Nx temperature constraints at points c^ 

(k = 1, 2, . . Nx), Nj Lagrangian multipliers are required, and equa- 
tion (22) of reference 3 can be rewritten for bar elements as 


- 


( 13 ) 


and for plate elements as 


Pn 


where 


Nt 

k=1 


Xk{'^}iCc]i{TQ}i^lj 

AiAi 


and 


"J =Z 


k=1 


Sit 




The resizing formulas become 


Ei\Y 

Ai,n+1 = [ — ) Ai,n 


for bar elements and 


'^j.n+l " i “ / 


for plate elements. 


(14) 


(15) 


( 16 ) 


(17) 


( 18 ) 
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In the resizing formulas the Lagrangian multipliers are determined by 
equation (7) of reference 4 rewritten as 




( 19 ) 


where Tjj and Tq are the temperature and the desired temperature at the 
kth constraint point, respectively, and b is a convergence parameter. Equa- 
tion (19) also permits consideration of inequality constraints of the form 


'^c,k 


( 20 ) 


since for passive constraints successive values of Xjj become small and eventu- 
ally approach zero and their impact in the resizing formulas becomes negligible. 

The number of iterations required to obtain a minimum-mass design depends 
on the value of the step-size parameter Y in equations (11), (12), (17), and 
(18) and on the value of the convergence parameter b in equation 09). Gener- 
ally, a trial -and-error technique is used to find values of Y and b which 
lead to minimum-mass designs with the fewest iterations. Since these parameters 
may be problem dependent, new values of Y and b should be determined for 
each new type of problem. To start the resizing process, an initial vector of 
design variables and initial values of are required in the resizing for- 

mulas. As a general procedure, for the problems considered herein, a unit 
initial design vector and unit initial values of are used to start the 

process . 


RESULTS AND DISCUSSION 

The effectiveness of the optimality criterion for problems with temperature 
constraints is assessed by applying the resizing formulas presented in the pre- 
vious section to the following problems involving uniform heating: 

(1) Cooling fin: A one-dimensional heat conduction problem with a single 

temperature constraint. Results for no convective heat loss are presented to 
demonstrate the effect of both the step-size parameter Y and the convergence 
parameter b on the resizing process and to assess effects of discretization 

on convergence to the closed-form solution for the minimum-mass design of cooling 
fins presented in reference 8. Effects of convective heat loss on the minimum 
fin mass for a given heating condition are also presented. 

(2) Square plate: A two-dimensional heat conduction problem with a single 

temperature constraint. Effects of convective heat loss on minimum plate mass 
for given heating conditions are considered, and results from the optimality 
criterion are compared with similar results obtained from the mathematical pro- 
graming technique of usable feasible directions described in references 9 and 10. 
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(3) Triangular plate; Similar to the square-plate problem except that 
convective heat loss is not included and effects of multiple temperature con- 
straints are considered. 

(4) Combined plate and bars: A problem with both one- and two-dimensional 

heat conduction paths and with multiple temperature constraints. Convective 
heat losses are not considered; however, two materials with different mass 
densities and thermal conductivities are used for the plates and bars to deter- 
mine a minimum-mass solution for a problem with mixed materials. 


Cooling Fin 

The geometry of a cooling fin with the associated boundary conditions is 
shown in figure 1 . Heat is generated in the fin at a constant rate of q per 
unit length, and convective heat loss may occur over both the top and bottom 
fin surface. The end X = 0 is maintained at a nondimensional temperature of 
zero, and the insulated end (point c) is constrained to a given value. The 
objective is to minimize the fin mass (or volume for constant material density) 
while satisfying the constraint condition. 

One -dimensional thermal finite elements with uniform areas and cubic tem- 
perature distributions were used to model the cooling fin to determine the tem- 
perature distribution. The elements have two nodes with T and dT/dx as 
nodal degrees of freedom. 

Results without convective heat loss .- The cooling fin problem without 
convective heat loss was used to determine values for the step-size parameter y 
and the convergence parameter b in the resizing formulas (eqs. (11) and (19)). 
Based on results from references 3 and 4, a value of 0.5 was initially used for 
y. The fin was idealized by five elements and minimum-volume solutions were 
£btained for a heat-generation number Nq = 1-0, a constraint temperature _ 

Tq = 0.3, and_values of b from 0.1 to 0.9. Results for the fin volume V, 
temperature Tq at the constraint point, and the number of iterations required 
to reach convergence are given in table I as functions of b. The number of 
iterations for convergence are also plotted as a function of b in figure 2. 

The data indicate that each solution converges to the same minimum fin volume; 
however, values of b from 0.5 to 0.7 require the fewest iterations to reach 
converged values of the fin volume and simultaneously satisfy the temperature 
constraint to four significant digits. The convergence parameter b = 0.5 
was selected as a best value for this problem. Additional solutions were then 
obtained for b = 0.5 and values of the step-size parameter y of 1 and 1/3. 
Results from these solutions revealed that y = 1 caused the resizing process 
to diverge and that y u "\/2> increased the number of iterations required for 
convergence. Since the other problems considered in this investigation were 
similar to the fin problem, values of y “ ^.5 and b = 0.5 were used for the 
remaining calculations. Values of these parameters, which produce most rapid 
convergence, however, are problem dependent and should be redetermined for 
problems which differ significantly from those considered herein. 


To study the effects of discretization on the finite-element procedure 
used for this problem, the fin was modeled with 5, 10, and 20 elements and 
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minimum fin’volumes were determined for Nq = 1.0 and Tq = 0.3. The results 
are given in table II and figure 3 and are compared with similar results from 
the closed-form solution of reference 8. Converged values of the minimum fin 
volume as a function of the number of elements are shown in figure 3(a). The 
finite-element results approach the closed-form solution slowly as the number 
of elements is increased; the difference between the two solutions decreases 
from 2.4_to 1.4 percent as the number of elements increases from 5 to 20. Fin 
volume V and constraint temperature Tq obtained from the 20-element solution 
are shown in figures 3(b) and 3(c), respectively, as a_function of iteration 
number. The data indicate that good convergence for V and Tq is obtained 
in 5 to 10 iterations. More detailed results are given in table II. 

Fin cross-sectional area distributions for the 5-, 10-, and 20-element 
solutions are compared with the area distribution from the closed-form solution 
in figure 4. The area distributions from the optimality criterion method slowly 
approach the closed-form solution as the number of elements increase. The 
20-element solution for the fin volume is within 2 percent of the closed-form 
solution. However, near the fin root, the areas given by the optimality crite- 
rion method are less than the closed-form solution. Near the tip, the areas 
given by the optimality criterion method are greater than the closed-form solu- 
tion. A finer mesh size near the fin tip or use of tapered elements could 
improve the agreement between the finite-element solution and the closed-form 
solution; however, since the objective of this investigation is to demonstrate 
the use of the optimality criterion, such refinements are considered 
unwarranted. 

The fin temperature distribution from the 20-element solution is compared 
with the closed-form solution in figure 5. At the center of the fin, the 
temperatures given by the finite-element solution are approximately 6 percent 
greater than the temperatures given by the closed-form solution. Again, this 
difference is attributable to use of constant-area elements and could be removed 
by use of a finer mesh size or tapered elements. 

Results with convective heat loss .- To determine the influence of terms 
in the heat-transfer matrix [D] , which are not functions of the design vari- 
ables (see discussion in the appendix), solutions for minimum fin volume were 
obtained by_using 5, 10, and 20 elements for values of from 0 to 0.5 and 

Nq = 1.0, Tq = 0.3, and T^ = 0.05. Converged values of minimum fin volume 
were obtained within 20 iterations for each idealization and value of B^. 

(See table III.) Since converged solutions were obtained within 20 iterations, 
effects on the resizing process of including constant terms in the heat-transfer 
matrix [d] appear to be negligible for this problem. Area and temperature 
distributions from the 20-element solution for B^ = 0 and 0.5 are shown in 
figures 6(a) and 6(b), respectively. The variation of minimum fin volume V 
with B^ is shown in figure 7. These results indicate that the minimum fin 
volume decreases and temperatures in the fin increase as convective heat losses 
increase. 
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Square Plate 


The effectiveness of the optimality criterion resizing formulas for two- 
dimensional problems with and without convective heat losses was assessed by 
obtaining minimum-mass (minimum volume for constant-mass density) solutions 
for a uniformly heated square plate with a single temperature constraint. Two- 
dimensional triangular thermal finite elements of constant thickness with a 
quadratic temperature distribution were used to determine the temperature 
response of the plate. The elements had six nodes, three at the vertices and 
three at the midpoints of the sides; temperature was the nodal degree of free- 
dom. The geometry, finite-element subdivision, and thermal boundary conditions 
for the plate are shown in figure 8. Temperature at point c was constrained 
to Tq = 10.0 and convective heat loss to the ambient atmosphere could occur 
from both top and bottom plate surfaces. Minimum-volume solutions were obtained 
for Nq = 1.0, Tg = 5.0, and values of = 0 (no convective heat loss), 
0.025, and 0.075. Iteration histories for the plate volume and temperature at 
the constraint point are given in table IV. Results from the mathematical pro- 
graming technique of references 9 and 10 are also shown. 

Results without convective heat loss .- Results in table IV indicate that 
in the absence of convective heat loss converged values of the plate volume and 
temperature at the constraint point are obtained with 20 iterations and agree 
with the mathematical programing results within approximately 1 percent. The 
thickness distribution corresponding to the minimum-volume solution for the 
finite-element mesh shown in figure 8 is presented in figure 9. Since the 
plate is symmetrical about the diagonal, only half of the plate is shown. The 
greatest material thickness is required along the insulated edges of the plate 
where the thickness increases from the constraint point to the heat sink. The 
temperature distribution in the plate for the minimum-volume solution is shown 
in figure 10 and indicates that, although the temperature constraint at point c 
is satisfied, the temperature at one point on the diagonal is approximately 

2- times the constraint temperature. Thus, if it is important that temperatures 
2 

over the entire plate do not exceed the value at point c, then temperatures 
at additional points must be constrained. Thickness and temperature distri- 
butions from the mathematical programing technique were virtually the same as 
those from the optimality criterion. 

Results with convective heat loss .- Iteration histories of V and T 
at the constraint point for the square plate for ® ®1 = 0.075 are 

shown in figures 11(a) and 11(b), respectively. In contrast to the cooling fin 
problem, where convergence occurred monotonically , convergence for the two- 
dimensional plate problem is oscillatory. Furthermore, the presence of the 
design-variable independent convective terms in the heat-transfer matrix causes 
the temperature response in the plate to become less sensitive to changes in 
the design variables and increases the number of iterations required for con- 
vergence. Results ^n table IV indicate that the optimality criterion yields 
values of V and T at the constraint point for B^ s 0.025 and 0.075 that 
agree within approximately 1 percent with similar values from the mathematical 
programing technique. 
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Triangular Plate 

The geometry, boundary conditions, and finite-element idealization for a 
uniformly heated triangular plate are shown in figure 12. The triangular ele- 
ments described previously were used to determine temperature response for the 
plate. Convective heat losses were not included in this problem. Minimum- 
volume solutions were obtained for the plate for Nq = 1.0 with temperature 
constraints at points c^ , C2, and C3. Two constraint combinations were con- 

sidered; namely, (1) point C3 wa^s constrained to Tq 3 = 10.0 with points c-| 
and C2 unconstrained and (2 ) Tq 3 = 10.0 with points c-| and C2 con- _ 
strained to Tq-| = Tq 2 = 5.0. Iteration histories for minimum plate volume V 
and temperatures at the constraint points from the optimality criterion for 
the two constraint combinations are shown in tables V(a) and V(b), respectively. 
Final results from the mathematical programing technique of reference 9 are also 
shown. A comparison of results from the optimality criterion and the mathemati- 
cal programing technique indicates that the two solutions for the minimum plate 
volume and constraint temperatures agree within about 1 percent. The iteration 
histories for plate volumes V and temperatures at the constraint points are 
shown in figures 13(a) and 13(b), respectively. Inclusion of the additional 
temperature constraints has two significant effects: (1) Since the constraint 

temperature at points ci and C2 is less than the corresponding temperature 
from the single constraint solution, additional plate volume is required to 
satisfy the constraints; and (2) the number of iterations required to reach a 
converged plate volume and simultaneously satisfy the multiple temperature con- 
straints is almost double that for the single constraint case. 

Temperature distributions for the single and multiple temperature con- 
straint conditions are shown in figures 14(a) and 14(b), respectively. Just 
as in the case of the square plate, when the temperature is constrained only 
at point C3, the temperature at other points on the plate greatly exceed the 
constraint temperature. The additional constraints at points c-| and C2 
improve but do not eliminate this condition. If it is important that temper- 
atures over the entire plate not exceed the value at point C3, temperatures 
at additional points must be constrained. 


Combined Plate-Bar Problem 

To illustrate further the usefulness of the optimality criterion, it was 
applied to the combined plate-bar problem shown in figure 15(a). Titanium was 
used for the plate material and aluminum was used for the bars. This structure 
was representative of a thermal design problem where one material satisfied 
strength requirements and the other acted as a highly efficient conductor to 
transfer the incident heat to a heat sink. The structure was uniformly heated 
at a rate of q = 34.1 kW/m^ (10 8OO Btu/ft^-hr) and was idealized as shown in 
figure 15(b) by eight of the previously described triangular thermal finite 
elements and four bar elements. Three noded constant-area bar elements with 
quadratic temperature distributions were used to model the bars to insure 
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compatibility between temperatures in the two types of elements. Temperatures 
in the structure were constrained not to exceed 260^ C (500^ F). Thus, based 
on the finite-element idealization, this problem had 20 temperature constraints. 

A converged solution for the minimum-mass configuration that satisfied the 
temperature constraints was obtained with 90 iterations. The large number of 
iterations required for convergence results from the large number of constraints 
in the problem. For the idealization used in the problem j a minimum mass of 
0.17 kg (0.37 Ibm) was obtained. When only the titanium plate was used, a 
minimum mass of 0.60 kg (1.33 Ibm) was required to satisfy the multiple temper- 
ature constraints. Thus, addition of the aluminum bars for the purpose of heat 
conduction reduced the structural mass by about 70 percent. 

The plate thickness distribution and bar area distribution for square bars 
corresponding to the minimum-mass configuration are shown in figure 16. Since 
the problem is symmetric, only half of the structure is shown. Nodal tempera- 
tures for the minimum-mass solution are given in table VI. A complete picture 
of the temperature field in the structure is shown in figure 17. Since the 
nodal temperatures presented in table VI are at discrete points, the tempera- 
tures shown in figure 17 are calculated over the entire plate at very small 
intervals. Again, because of the symmetry in the problem, temperatures are shown 
for only half of the structure. The sharp dips in the temperature distribution 
near the aluminum bars occur because the bars provide a better conduction path 
to the heat sink than the titanium plate. 


CONCLUDING REMARKS 

An optimality criterion is presented for minimum-mass design of heated 
structures with constraints on temperature. The thermal optimality criterion 
is similar to an existing optimality criterion for mechanically loaded struc- 
tures with displacement constraints- Resizing formulas from the thermal opti- 
mality criterion are applied to both one- and two-dimensional problems where 
temperatures can be controlled by varying material distribution in the structure. 
A comparison of results obtained from the optimality criterion with results from 
a closed-form solution and with results from a mathematical programing technique 
reveal that the optimality criterion results agreed within 2 percent of those 
from the other methods. 

For the one-dimensional problems the optimality criterion converges mono- 
tonically. However, for the two-dimensional problems convergence occurs in 
an oscillatory manner, and the number of iterations required for convergence 
increases significantly with either the number of temperature constraints or 
degree of convective heat loss present in the problem. 

The thermal optimality criterion is used to size components of a heated 
structure composed of titanium plates and aluminum bars with temperatures con- 
strained to 260^ C (500^ F) or less. Based on thermal considerations only, 
adding aluminum bars reduces the combined mass by about 70 percent of that for 
titanium plates alone. 
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The thermal optimality criterion augments optimality criteria for mechan- 
ically loaded structures with strength and stiffness related constraints and 
offers the possibility of extension of optimality techniques to resizing struc- 
tures with combined thermal and mechanical loading. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
August 9 , 1 977 


APPENDIX 


THERMAL OPTIMALITY CRITERION 
Single Temperature Constraint 

For structures with both conductive and convective modes of heat transfer 
present, the equation governing temperature response in the structure can be 
written in finite-element matrix notation as 


[D] {t} = {q}i + {q}2 


(A1) 


where [I^ is a heat-transfer matrix composed of terms for both conductive and 
convective heat transfer, {t} is a vector of unknown temperatures, and {q}-| 
and {q }2 are applied thermal and convective heat load vectors, respectively. 
The matrix [dJ can be written as 

[d] = [Cj + [h] (A2) 

where [c] is a conduction matrix and [H] is a convection matrix. For struc- 
tures which can be represented by N 5 bar elements of constant area and Np 
plate elements of constant thickness, the total mass can be expressed as 

Nb Np 

i =1 j =1 


where P 5 and Pp are the mass densities of the bar and plate materials, 
respectively; and are the cross-sectional area and length of the 
ith bar element, respectively; tj and Sj are the thickness and surface 
area of the jth plate element, respectively. The design variables are A^ 
and t j . 

The temperature at some point c on the structure is constrained to a 
specified value Tq. The equation of constraint is 



m 


I 


Pb^i^i 


T - Tq = 0 


(A4) 


or 


(A5) 


{t}T{q} - Tc = 0 
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where \T} is the nodal temperature vector obtained by solution of equation (A1) 
and {q} is a vector of zeros except for a value of unity at the node corre- 
sponding to point c. The optimality criterion is given by equation (4) in the 
main text and is rewritten in terms of the thermal problem as 


^ -T— fo} = 0 

9Ai 


(A6) 


for bar elements and 


PpSj * X ° 


(A7) 


for plate elements where X is an unknown Lagrangian multiplier. 
The vector {q} can be written as 


{q} = [d]{Tq} (A8) 

where {Tq} is a temperature vector obtained by solution of equation (A1) for 
a unit thermal load applied at point c. Replacing {q} in equations (A6) and 
(A7) with its equivalent from equation (A8) gives 


9 { t} 

Pb^i + ^ {tq} = 0 

o Aj^ 


and 


3{t}T 

PpSj + ^ -^[d]{tq} = 0 


3t^ 


(A9) 


(A10) 


as expressions of the optimality criterion for bars and plates, respectively. 

Expressions for the partial derivatives of the temperature vector can be 
obtained by differentiation of equation (A1). Note that, for the type of ther- 
mal finite elements used to model the structure, only the conduction matrix [c] 
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is proportional to and Differentiation of equation (A1) with respect 

to A^ gives 


3Ai ^ ^ 3Ai 


0 


(All) 


or 


9{t}T 

3Ai 





9[cJi 

3Ai 


(A12) 


Similarly for tj 


3{t}T 

It” 


[»] 



3[C]j 

3 1 j 


(A13) 


Since, for the type elements used herein, the conduction matrix [cfj is pro- 
portional to Aj^ and tj, the partials of [Cj i with respect to A^ and tj 
can be written as 


3[c]i [Cli 

3Ai ■ ~ 


(A14) 


and 


3 [c] j j 


(A15) 


Thus, the optimality criterion can finally be written as 


T i 

Pbi^i - X{T}i — {Tq}i = 0 
*i 


(A16) 
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for bars and 


r ,T [C]j. 

PpSj - X{T]j — {Tq}j = 0 


(A17) 


for plates. The subscripts i and j on the quantities {t}, [cJ , and {Tq} 
indicate that those quantities are associated with the ith bar element or jth 
plate element. 


Multiple Temperature Constraints 

If the structure has Nj temperature constraints at points cj^ 

(k = 1, 2, . . N't), the constraint equations are written as 

{T}T{q}k - To^k =0 (k = 1, 2 Nt) (A18) 


The functional given by equation (3) in the main text becomes 


Nb Np Nt 

m» PpAiJli + '2 PptjSj + 2 ^k[{T}T{q}k - Te,k] 

i=1 >1 k=1 


(A19) 


where are the unknown Lagrangian multipliers for each constraint k. 

Application of the same procedure used for a single constraint leads to the 
optimality criterion for multiple constraints 


Np 


[C]i, 


Pb% “ 2 ^ ^k't'^^i i ,k = ® 


(A20) 


for bars and 


Nn 


T 0^3 j 

k=1 ^ 


(A21) 


for plates. 
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TABLE I.- EFFECT OF b ON CONVERGENCE OF RESIZING PROCESS 
FOR COOLING FIN PROBLEM IDEALIZED WITH FIVE ELEMENTS 

[Nq = 1.0; Bi = 0; % = O.3] 


b 

Converged 
value of volume 
V 

Converged value 
of temperature 
at constant point 
T 

No. of iterations 
to reach convergence 

0.1 

1 .5168 

0.3000 

139 

.3 

1 .5168 

.3000 

38 

.5 

1 .5168 

.3000 

16 

.7 

1 .5168 

.3000 

16 

.9 

1 .5168 

.3000 

23 


TABLE II.- EFFECT OF FINITE-ELEMENT MESH ON CONVERGENCE OF 
RESIZING PROCESS FOR COOLING FIN PROBLEM 
|Ng = 1.0; Bi = 0; Tq = 0.3; Y = 0.5; b = 0.^ 


Iteration 

5 elements 

10 elements 

20 elements 

V 

Tc 

V 

Tc 

V 

Tc 

1 

0.6719 

0.6691 

0.6686 

0.6674 

0.6674 

0.6669 

5 

1.4441 

.3148 

1 .4342 

.3146 

1 .4306 

.3146 

10 

1 .5145 

.3004 

1.5039 

.3004 

1 .5001 

.3004 

15 

1 .5167 

.3000 

1 .5061 

.3000 

1 .5022 

.3000 

20 

1.5168 

.3000 

1 .5061 

.3000 

1 .5023 

.3000 

(a) 

1 .4815 

.3000 

1 .4815 

.3000 

1 .4815 

.3000 


^Obtained from closed-form solution of reference 8 





TABLE III.- EFFECT OF CONVECTION ON MINIMUM VOLUME 
FOR COOLING FIN 

[Nq = 1.0; Tq = 0.3; Tg = 0.05; 20 iterationsj 


No. of 
elements 

Converged fin volume, V 

Bi = 0 

Bi = 0.05 

Bi = 0.2 

Bi = 0.5 

5 

10 

20 

1.5168 
1 .5061 
1 .5023 

1 .4878 
1 .4774 
1 .4736 

1 .4005 
1 .3906 

1 .3870 

1 .2232 
1 .2147 
1 .2115 
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TABLE IV.- EFFECT OF CONVECTION ON MINIMUM VOLUME FOR HEATED 


SQUARE PLATE WITH SINGLE TEMPERATURE CONSTRAINT 
[ng = 1.0; = 10.0; fg = 5.o] 




Bi = 0 



□ 

).025 


Bi = 

0.075 


Iteration 

Plate 

Temperature 









volume 

at constraint 


V 


T 

V 


T 




point 








0 

1 .0000 

0 

.59 

1 .0000 

0 

.66 

1 .0000 

0 

.78 

1 

.7521 


.74 

.7890 


.78 

.8530 


.86 

2 

.1787 

2 

.99 

.1970 

2 

.90 

.2330 

2 

.78 

3 

.0478 

10 

.70 

.0523 

9 

.42 

.0618 

7 

.77 

4 

.0259 

19 

.80 

.0246 

16 

.70 

.0242 

12 

.70 

5 

.0268 

18 

.80 

.0199 

19 

.00 

.0135 

15 

.20 

6 

.0375 

13 

.30 

.0239 

16 

.50 

.0097 

16 

. 10 

7 

.0512 

9 

.65 

.0347 

12 

.30 

,0091 

16 

.10 

8 

.0588 

8 

.31 

.0491 

9 

.22 

.0115 

15 

.10 

9 

.0574 

8 

.41 

.0581 

7 

.91 

,0134 

13.00 

10 

.0521 

9 

.18 

.0571 

7 

.93 

.0324 

10 

.20 

11 

.0475 

9 

.95 

.0505 

8 

.76 

.0510 

7 

.93 

12 

.0455 

10 

.30 

.0442 

9 

.80 

.0625 

6 

.95 

13 

.0453 

10 

.30 

.0405 

10 

.50 

.0598 

7 

.15 

14 

,0461 

10 

.20 

.0395 

10 

.70 

.0492 

8 

.08 

15 

.0470 

10 

.00 

.0403 

10 

.60 

.0388 

9.28 

16 

.0476 

9 

.91 

.0419 

10 

.20 

.0317 

10 

.30 

17 

.0477 

9 

.91 

.0434 

9 

.90 

.0278 

11 

.00 

18 

.0476 

9 

.95 

.0443 

9 

.80 

.0263 

11 

.30 

19 

.0474 

10 

.00 

.0444 

9 

.89 

.0269 

11 

.20 

20 

.0474 



.0444 

10 

.00 

.0289 

10 

.30 

21 

.0473 



.0444 

10 

.00 

.0318 

10 

.30 

22 

.0474 



.0444 

10 

.00 

.0347 

9 

.89 

23 





.0444 

9 

.98 

.0368 

9 

.58 

24 





.0442 

9 

.95 

.0378 

9 

.46 

25 





.0439 

9 

.95 

.0374 

9 

.53 

26 





.0436 

9 

.97 

.0363 

9 

.71 

27 





.0434 

10 

.00 

.0350 

9 

.92 

28 





.0434 



.0339 

10 

.10 

29 





.0435 



.0332 

10 

.20 

30 





.0437 



.0330 

10 

.20 

31 





.0439 

\ 


.0333 

10 

.20 

32 





.0440 

9 

.98 

.0337 

10 

.10 

33 





.0439 

9 

.97 

.0342 

10.00 

34 





.0438 

9 

.97 

.0346 

9 

.95 

35 





.0437 

9 

.99 

.0348 

9 

.91 

36 





.0436 

10 

.00 

.0349 

9 

.91 

37 





.0436 

j 


.0347 

9 

.93 

38 





.0436 



.0345 

9 

.90 

39 





.0437 



.0343 

10 

.00 

40 





.0438 

1 

j 


.0341 



41 





.0438 

i 


.0340 



42 





.0438 

9 

.99 

.0340 



43 





.0438 

9 

.99 

.0341 



44 





.0437 

10 

.00 

.0342 



45 









.0343 



46 









.0344 

9 

.99 

47 









.0344 

9 

.98 

48 









.0344 

9 

.98 

49 









.0344 

9 

99 

50 









.0343 

10 

.00 

51 

1 







f 

.0343 

10 

.00 

(a) 

.0479 

9 

.98 

.0442 

9 

.98 

.0342 

9 

97 


^Obtained by method of references 9 and 10. 



TABLE V.- EFFECT OF MULTIPLE TEMPERATURE CONSTRAINTS 
ON HEATED TRIANGULAR PLATE 
[Nq = 1.0; Bi = O] 


(a) Single constraint 

[to3 = 'f-O] 


Iteration 

V 

Tel = Tc2 

Tc3 

0 

0.5000 

0.38 

0.53 

1 

.5140 

.36 

.51 

2 

.1180 

1 .57 

2.21 

3 

.0266 

6.97 

9.76 

4 

.0125 

14.80 

20.70 

5 

.0123 

15.00 

20.90 

6 

.0176 

10.50 

14.50 

7 

.0255 

7.25 

10.00 

8 

.0307 

6.03 

8.31 

9 

.0307 

6.03 

8.29 

10 

.0279 

6.63 

9.08 

11 

.0254 

7.29 

9.97 

12 

.0242 

7.66 

10.40 

13 

.0241 

7.69 

10.40 

14 

.0247 

7.54 

10.20 

15 

.0252 

7.38 

9.99 

16 

.0255 

7.32 

9.87 

17 

.0254 

7.34 

9.87 

18 

.0252 

7.40 

9.93 

19 

.0251 

7.47 

9.99 

20 

.0250 

7.41 

10.00 

21 

.0249 

7.41 

10.00 

22 

.0249 

7.45 

10.00 

23 

.0250 

7.47 

9.99 

24 

.0249 

7.43 

10.00 

25 

.0249 

7.43 

10.00 

(a) 

.0252 

7.41 

9.99 


^Obtained by method of references 9 and 10. 




TABLE V.- Concluded 


(b) Multiple constraints 

[fc 3 = 10.0; = fc 2 = 5.o] 



^Obtained by method of references 9 and 10. 


2U 


















(a) 5-element solution. 



(b) 10-element solution. 



(c) 20-element solution. 

Figure 4.- Area distribution_for fin problem. 
Nq = 1.0; Bi D 0; Tq = 0.3. 
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00 


Figure 10.- Temperature distribution for minimum volume square plate. 





(a) Convergence of volume. 



(b) Convergence of temperature at constraint point. 

Figure 11.- Effect of convection_on square plate. 
Nq = 1.0; Tq = 10.0; Ta = 5.0. 
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.04 

.02 


0 10 20 30 40 

No, of iterations 

(a) Convergence of volume. 



0 10 20 30 40 

No. of iterations 

(b) Convergence of temperature at constraint points. 

Figure 13.- Effects of multiple constraints on triangular plate 

Nq = 1.0; Bi = 0. 





(b) Multiple constraints. 

Figure 14.- Temperature distributions for heated triangular plate. 
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(a) Geometry and material properties. 



(b) Finite-element idealization. 

Figure 15.- Combined plate-bar problem. Heat flux, 3^.1 kW/m^ 
(10 800 Btu/ft^-hr) . (Dimensions in cm (in.).) 
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0,91 (0.36) 



Figure 16.- Material distribution for combined plate-bar problem. 

(Dimensions in cm (in.).) 
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